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ABSTRACT 



It has long been realized that dark matter halos formed in cosmological N-body 
simulations are characterized by density profiles p{r) that, when suitably scaled, have 
similar shapes. Additionally, combining the density and velocity dispersion profiles 
^ ' '^C'^)) each of which have decidedly nonpower-law shapes, leads to quantity p/a^ that is 

a power-law in radius over 3 orders of magnitude in radius. Halos' velocity anisotropy 
profiles /3(r) vary from isotropic near the centers of halos to quite radially anisotropic 
near the virial radius. Finally, there appears to be a nearly linear correlation between 
(3 and the logarithmic density slope 7 for a wide variety of halos. This work is part 
of a continuing investigation of the above interrelationships and their origins using 
analytical and semi-analytical techniques. Our finding suggest that the nearly linear 
/3-7 relationship is not just another expression of scale- free /o/o"^ behavior. We also 
note that simultaneously reproducing density and anisotropy profiles like those found 
in simulations requires /3(r) and 7(r) to have similar shapes, leading to nearly linear 
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/9-7 correlations. This work suggests that the j3-^ and power-law p/a^ relations have 
distinct physical origins. 

Subject headings: dark matter — galaxies: structure — galaxies:kinematics and dynam- 
ics 

1. Introduction 

Cosmological N-body simulations are important tools that allow us to reconstruct how gravi- 
tationally bound objects like galaxies and clusters form and evolve in the universe. The base-level 
simulations involve Newtonian gravitational interactions between particles that represent masses of 
dark matter exclusively. These types of simulations suppress the short-range interactions between 
masses, making the dynamics collisionless, i.e., particles move under the influence of a global po- 
tential. The consensus view of the results of such simulations is that all bound structures have 
common density structures determined by a few parameters. The most often discussed density 
distribution in the literature is the empircal Navarro-Prenk- White (Navarro, Frenk, &: White 1996, 
1997, NFW) function. However, many different functions have been invoked to describe the density 
profiles of dark matter structures resulting from individual simulations {e.g., Sersic 1968; Moore 
et al. 1998; Navarro et al. 2004). While there is a considerable amount of literature attempting to 
explain the universality of N-body density profiles (Syer & White 1998; Lokas 2000; Nusser 2001; 
Barnes et al. 2005; Lu et al. 2006), no widely accepted explanation exists. 

We are looking to gain a fundamental understanding of the physics driving the universality 
of dark matter halos in collisionless simulations. Such insight may eventually provide a way of 
reconciling the differences between simulated halos and the inferred properties of halos of observed 
galaxies, such as the cusp-core controversy {e.g., de Blok, McGaugh, & Rubin 2001; Spekkens, 
Giovanelli, k. Haynes 2005). However, the results of the present paper cannot be directly applied 
to observed galaxies and clusters because the latter are affected by dissipational baryonic pro- 
cesses, which we ignore in the present work; we concentrate exclusively on the effects of gravity in 
collisionless systems. 

In the past few years, two interesting relationships have been uncovered that link various 
dynamical properties of halos formed in simulations like those discussed above. Taylor &: Navarro 
(2001) have found that there is a scale-free relationship between a phase-space proxy and radius; 
p/cT'^ oc r~", where p is density, a is velocity dispersion, r is radial distance, and a is a constant for 
any given halo. Interestingly, this relationship is not unique to N-body results. Austin et al. (2005) 
have shown that halos created semi-analytically also have power-law p/cr^ behavior and speculate 
that this feature is a result of violent relaxation. 

The other relationship relates the density and velocity anisotropy profiles of N-body halos 
(Hansen & Moore 2004). Specifically, the curve describing the logarithmic radial derivative of the 
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logarithm of the density has the same basic shape as the radial profile of a halo's anisotropy. As 
with the power-law behavior of pjcr' ^ there is no obvious basis for this relationship, but together 
they signal that there are unifying physical processes at work in halo formation. 

This paper and its companion (Barnes et al. 2006, hereafter Paper I) take steps towards 
uncovering what these processes are. In that work, we found that the mechanical equilibrium of 
halos with velocity isotropy is not the cause of power-law p/a^ profiles, as at least two well-known 
equilibrium density profiles (Hernquist and King) do not produce p/cr^ power-laws. However, 
Navarro et al. (2004) and Sersic density distributions in mechanical equilibrium produce nearly 
scale- free p/cr^ profiles. We have also investigated density profiles that result from solving the 
isotropic Jeans equation with the constraint that p/cr^ oc r~". In this restricted case, density 
profiles are similar to those presented in Navarro et al. (2004) but are reproduced better by Sersic 
models. 

We continue and expand this work by again considering the density profiles that result from 
solving the Jeans equation under the assumption of scale- free p/cr'^ (motivated by N-body and 
semi-analytical results), but now the velocity distributions will not be restricted to be isotropic. 
The details of solving the "constrained" Jeans equation are given in §2. The anisotropic Jeans 
equation admits a wide variety of solutions. We use generic, but sensible, forms for our adopted 
anisotropy profiles. These profiles and the solutions that result are discussed in detail in §3 and 4. 
Our main conclusion follows from these sections; equilibrium density profiles replicate those seen 
in simulations if the radial density and anisotropy profiles of any given halo have similar shapes. 
The degree of connectedness between the p/a^ and density slope-anisotropy relations is considered 
in Section 5. We end with a summary of our findings and our conclusions. 



2. The Constrained Jeans Equation 

The Jeans equation is a moment of the collisionless Boltzmann equation. Where the Boltz- 
mann equation describes the evolution of the full phase-space distribution function, which is not 
easily obtained from simulations, let alone observations of real galaxies, the Jeans equation relates 
observationally accessible quantities, like density and velocity dispersion. Mechanical equilibrium 
for a spherical collisionless system with radially dependent anisotropy is determined through the 
Jeans equation (Jeans 1919; Binney & Tremaine 1987), 



d 
dr 



p{r)a'^{r) 
3 - 2/3(r) 



, 2(5{r) p(r)a2(r) _ „ , .M(r) 



where M(r) is the mass enclosed at radius r, p is the density, a is the total velocity dispersion, and 
j3 is the anisotropy. The 3 — 2(5{r) terms result from our choice to use the total velocity dispersion 
rather than the r-component. While we do not directly calculate phase-space distribution functions 
here, the density and velocity dispersion functions that will be discussed are everywhere positive, 
implying that the corresponding distribution functions are physically plausible. 
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One can reduce the number of functions in Equation 1 by assuming an "equation of state" 
that connects the density to the total dispersion. We use the empirically established relation 
pja^ = {po/vQ){r/rQ)~°'. We note that while this connection is consistent with the available 
numerical evidence (Taylor & Navarro 2001; Dehnen & McLaughlin 2005), it is not the only choice. 
One could equally well choose to utilize the radial dispersion, as in Dehnen &; McLaughlin (2005). 
At this point in time, it is not clear which choice is the most physically relevant. 

Imposing this scale- free p/cr^ constraint and changing to the dimensionless variables x = r/r^ 
and y = p/ po, we rewrite Equation 1 as. 



x^ 






y 


dx 


l^3-2/3(x) ) 



+ 



- W{x) 



y5/3^2a/3- 



BM{x), 



where B = G/tqVq. Differentiating this equation with respect to x gives us. 



d_ 

dx 




y5/3^2a/3-l 



Cyx\ 



(2) 



(3) 



where C = 47rpo'"o/^o- This expression is an extension of equation (8) in Paper L We eliminate the 
constant C by solving for y, differentiating with respect to x again, and grouping like terms. The 
resulting constrained Jeans equation is. 



(2a + 7 - 6)(-(q - 7) + l)(2a - 57) 



157" + 37'(8a - 57 + 4/3 + 12/356i - 5) 
-3(5[6i(4a2 + 7^ - 807 + 8q + 77 - 15)] 



(4) 



-3<52 [661 62(0 -7 + 1)] 
-3(53[63(54/3 + 144/3^ + 24/?^)] 

-3(5' [661 (a - 7 + 1) + 961625] - 3(5"[36i] 

In this notation, 7 = 7(x) = —d\ny/d\nx is the logarithmic density slope, 6 = 5{x) = dln/J/dln x, 
61 = 2/3/(3-2/?), 62 = (3 + 2/3)7(3-2/3), 63 = (3-2/?)-3^ and the primes indicate derivatives with 
respect to Inx. 

As discussed previously, the parameter a is a constant, and we choose to investigate the 
following values; a = 1.875, 35/18 = 1.944, 1.975. These values have been chosen based on the 
findings of Austin et al. (2005). For isotropic systems, a = 35/18 represents a critical point for 
solutions of the constrained Jeans equation. Solutions with a > 35/18 have 7 values that asymptote 
to finite values, and a = 1.975 is our example of such a case. In terms of density, these solutions 
approach power-law density profiles at large radius. In the other case, 7 values increase indefinitely 
with radius, implying density profiles that continually steepen farther from the center. Our value 
of a = 1.875 (which is also the value found by Taylor &: Navarro (2001)) represents this class of 
solutions. 



While aesthetically unappealing in the extreme, Equation 4 can be straightforwardly solved 
numerically for 7 over a range of Inx once a and /3(x) have been specified (different choices of P{x) 
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are considered in §§3, 4.1, and 4.2). Our approach to the solution sets 7 = 2 at the scalelength ro 
(x = 1). In effect, we define the scalelength of a halo to be the radius where the density slope has 
the isothermal value. We further define Xvir = 10, by analogy with c = 10 NFW profiles [note that 
this does not mean that p{x < 10) = 200pcrit] ■ The remaining initial condition is the value of 7' at 
tq. We leave this as a free parameter to be determined by minimizing the differences between the 
solution 7-profile and those of the NFW, 

pix)/po = x-^(l + x)-2 

7(x) = (l + 3x)/(l+x) (5) 

and Navarro et al. (2004, N04) 

log[p(x)/po] = (-2/^)(x'^ - 1) 

7(x) = 2x'' (6) 

density profiles. In particular, for the NFW profile 7'(1) = 0.5 and for the N04 profile 7'(1) = 2/i, 
where ^ = 0.17 (the best-fit value from N04). We have chosen these profiles because of their ability 
to empirically describe the collisionless halos that form in a wide variety of numerical simulations 
as well as their relative simplicity, but any other profile could be chosen and analyzed in the same 
way. 



3. Constant Anisotropy Solutions 

The simplest type of anisotropy distribution is a constant nonzero value. In this case, equation 4 
simplifies considerably as 5 and its derivatives disappear. Our work in Paper I deals with the 
special case of /? = 0. Here, we extend to nonzero, but constant, anisotropics (—0.2 < (3 < 1.0). 
The relevant constrained Jeans equation is, 

2 

(2a + 7 - 6) (-(a - 7) + l)(2a - 57) = I57" + 37'(8a -57 + 4/3-5). (7) 

Solutions to this equation extend and agree with our earlier work. In general, solutions to 
the constant anisotropy constrained Jeans equation do not resemble NFW profiles for any tested 
combination of constant anisotropy value and a value. On the other hand, the constrained Jeans 
equation solutions can approximate N04 profiles, at least for isotropic and slightly tangentially 
anisotropic velocity dispersion distributions. For the solutions that best approximate N04 profiles, 
the anisotropy value decreases as the a value increases. Figure 2 shows 7 profiles for a = 1.875, j3 = 
0.1 (panel a), a = 35/18 = 1.944, /3 = -0.1 (panel b), and a = 1.975, /3 = -0.2 (panel c). Even 
mild radial anisotropy {(3 > 0.3) drives the solutions to have fairly constant 7 profiles, independent 
of a. Panel d of Figure 2 shows the 7 profile for a = 1.875, (3 = 0.5. We note that the special value 
a = 35/18 is obtained analytically only for /? = (Austin et al. 2005). For nonzero constant P, the 
special value of a can be found by setting the 7' term in Equation 7 equal to zero and utilizing the 
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relationship for power-law density profiles, 7pi = 6 — 2a. The result is that agpcciai = 35/18 — 2/5/9 
(see also Dehnen &: McLaughlin 2005). 



4. Solutions for Anisotropy Distributions 

On relaxing the constant anisotropy requirement, we are immediately faced with a decision. 
What form should the anisotropy distribution P{x) take? We have chosen two flexible functions to 
investigate; one is a simple power-law while the other is based on the hyperbolic tangent function. 
For both functions, we choose boundary anisotropy values and we have parameters that affect 
the rates of change between these boundary values. The functions are similar in that they both 
monotonically increase with Inx, following the stability argument given in Merritt &; Aguilar (1985). 
As the next two subsections will show, the important difference between the two functions is that 
the power-law has a positive second derivative everywhere while the hyperbolic tangent has an 
inflection point. 

Specifically, the power-law function is given by, 

Inx + lnxbound 



/3(x) = Ao + {Phi - Ao) 



In Xvir ~l~ 111 27i-|ound 



(8) 



where Ao is the anisotropy value as In x — > — In Xbound , Phi is the anisotropy value at the scaled 
virial radius Xvir, ij > controls the shape of the distribution, and the range of integration is 
determined by Inxbound = 25. We fix the scaled virial radius at the value for an NFW halo with 
a concentration index of 10, x^ir = 10. While our solutions cover large ranges of x, we will be 
interested primarily in the several orders of magnitude surrounding the scalelength. 

The other distribution that we study will be referred to as the tanh anisotropy profile. The 
expression for this distribution is, 

Pix) = ^(/3hi-Ao)[l+tanh(slnx)]+Ao, (9) 

where again Ao is the anisotropy as Inx —> — Inxbound, Phi is now the anisotropy as Inx ^ Inxbound, 
and s determines the rate of change between these boundary values. Now, tanh z = (e^ — e~^) /(e^ + 
e~^) and 1 + tanhz = 2e^/(e^ + e~^). If z = slnx = Inx**, then 1 + tanhz = 2x**/(x* + x~*) = 
2x?^ I{x?^ + 1). The tanh form is thus a generalization of commonly discussed anisotropy functions; 
P (X r"^ /{r"^ + r^) (Osipkov-Merritt) for s = 1, and P cc r/{r + a) {e.g., Mamon & Lokas 2005) for 
s = 1/2. The tanh form we adopt allows for more flexibility in the anisotropy profile. 

Our approach to solving Equation 4 with these anisotropy distributions is as follows. First, 
choose a profile type to fit to; either the NFW or N04 7 distribution. Next, choose pairs of (Ao; 
Phi) values. With the boundary anisotropy values fixed, we are left with one free parameter in each 
P profile. For each set of (Aoj Phi)^ we find the r] 01 s that minimizes the average least-squares 
difference between the 7 profile that results from solving Equation 4 and the chosen 7 distribution 
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over the range 10~^3;vir and Xvir, a range over which N-body simulations produce meaningful results. 
This procedure provides us with the best approximation to either the NFW or N04 profiles given a 
set of boundary anisotropy values. Note that we do not use the least-squares values for hypothesis 
testing or for making statistical comparisons between the various anisotropy forms, but rather to 
find the form of the anisotropy profile that best reproduces NFW/N04 density profiles, while solving 
the Jeans equation. 

4.1. Power-law Anisotropy Solutions 

For the power-law anisotropy solutions, we have investigated a parameter space with 0.0 < 
Ao ^ 0.9 and 0.1 < /3hi < 1-0. Points where /?io = /3hi have not been included as they are 
identical to the constant anisotropy models (§3). Points with Aoi/^hi < have been excluded due 
to an instability in the fitting routine for these models. We have also investigated the regions of 
parameter space where /3io > /3hi, but those fits are uniformly uninteresting as the solutions do not 
resemble any sort of physically relevant density profile. The solutions which best match NFW and 
N04 profiles all have /3io = 0. Solutions that are matched to NFW profiles tend to have smaller 
best-fit /3hi values than those matched to N04 profiles. The best-fit solutions are discussed in detail 
below. 

The best-fit 7 distributions with power-law anisotropy profiles are shown in Figure 3, where the 
solutions (solid lines) are matched to the NFW profile (dashed lines). Dash-dotted lines illustrate 
N04 7 profiles. Dotted lines show the anisotropy profile from Mamon & Lokas (2005). Panels 
a and b show the 7 and /? profiles calculated with a = 1.875 that best match the appropriate 
profile, panels (c, d) and (e, f) illustrate the same distributions, but for a = 35/18 and a = 1.975, 
respectively. Note that the solutions for the power-law anisotropy profile never resemble the NFW 
profile; i.e., the best-fits are not good fits. 

Figure 4 illustrates the same relations, but the solutions are matched to the N04 profile. The 
line styles are the same as in Figure 3. The N04 profile is reproduced with much greater accuracy 
than the NFW profile, but the match becomes substantially worse as a increases. Note that for the 
solutions that match the N04 profiles, the best-fit anisotropy profile has an isotropic distribution 
near the center and a highly radial anisotropic character (0.8 ^ /9 ^ 1.0) near the virial radius. 
This character of anisotropy profiles has been noted previously in N-body simulations {e.g., van 
Albada 1982; Cole & Lacey 1996). It has also been suggested that such behavior is evidence that 
an orbital stability process is at work in determining the density structure of equilibrium halos 
(Barnes et al. 2005). 

In numerous N-body simulations, a nearly linear relationship has been found between the 
anisotropy and density-slope profiles (Hansen & Moore 2004, 2006). It is important to note that 
no one simulation displays an exactly linear /3-7 profile. Over the range of 7 values our solutions 
cover (1 ^ 7 ^ 3), individual simulation fd-j curves are nearly linear indicating that the anisotropy 
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and density slope distributions have similar shapes. The variations present in Figure 2 of Hansen 

6 Moore (2006) appear to originate from differing slopes and intercepts for different simulations. 

In Figures 5 and 6, we present the (3 vs. 7 profiles (solid lines) obtained from the solutions 
shown in Figures 3 and 4 and compare them to the relations found in N-body simulations. The 
dashed lines represent the mean trend of the simulation results and the dash-dotted lines show 
the extent of variations in those results. Of the six curves in Figures 5 and 6, the best qualitative 
match to the near-linear relationship found in Figure 2 of Hansen & Moore (2006) is given by the 
N04-like profile with a = 1.875, but other values 1.8 < a < 1.9 would also give similar results 
(Figure 6a). One difference is that our relation is "steeper"; we reach total radial anisotropy at 

7 = 3, while the mean simulation results reach (3 ~ 0.5 at 7 = 3. Solving the constrained Jeans 
equation for N04-like profiles with a > 1.9 or for NFW-like profiles does not produce near- linear 
/3-7 relations, but the solutions are at least consistent with the mean trend within the variations 
seen in simulations. 



4.2. Tanh Anisotropy Solutions 

We now turn to the hyperbolic tangent (tanh) anisotropy distribution solutions. The parameter 
space for the tanh profile is slightly different than for the power-law anisotropy distribution; —0.2 < 
Ao < 0.9 and 0.3 < /3hi < 1-0. Other combinations of Aoand /3hi yield unrealistic density profiles. 
As in the power-law case, the Ao values for the solutions that best match both NFW and N04 
profiles are small, Ao < 0.1. The /3hi values for the NFW matched solutions are larger than those 
for the N04 matched solutions. The specifics of the best-fit solutions are given below. 

The solutions that best approximate the NFW profile have anisotropy profiles that vary from 
isotropic near the center to mildly radially anisotropic {(3 > 0.5) at the virial radius. The best-fit 
solutions' 7 and (5 profiles are shown in Figure 7. Again, the line styles are the same as in Figure 3. 
Unlike in the power-law (3 case, it is easy to see that these solutions approximate the asymptotic 
NFW behavior quite well. 

Another difference from the power-law case is the behavior of the solutions that resemble 
N04 profiles. These solutions tend to have fairly constant anistropy distributions that are nearly 
isotropic {(3 < 0.3). This is most evident when we look at the best-fit solution 7 and (3 profiles 
in Figure 8. The N04 7 profile is well-fit by the solution, but only when the anisotropy profile is 
nearly constant near zero. This agrees with the work presented in Paper I where a /3 = profile 
produces solutions that are good approximations to N04 profiles. 

Figures 9 and 10 show the (3 vs. 7 curves for the NFW-like and N04-like solutions, respec- 
tively. The curve for the NFW-like solution with a = 1.875 (Figure 9a) is remarkably linear, again 
indicating the similarity of the (3 and 7 distributions. In fact, each of the curves in Figure 9 have 
substantial regions of linearity. This behavior makes them all qualitatively similar to the results 
of Hansen & Moore (2006) (marked by the dashed and dash-dotted lines as before). However, the 
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most linear relation (panel a with a = 1.875) is also the one that disagrees most strongly with 
the simulation results. The remaining profiles {a = 1.944 and a = 1.975) actually show a decent 
quantitative match to the mean simulation trend. Figure 10 shows that the N04-like density profile 
that is most linear (again a = 1.875) is a poor quantitative match to simulation results. Overall, 
the N04-like density profiles produce nearly horizontal (3-'y profiles, qualitatively dissimilar to the 
profiles from simulations. However, we note that profiles created with larger (> 1.944) a values are 
quantitatively consistent with simulations. 

5. Relating Scale- free p/a^ and f3 vs. 7 

As these results have been found by assuming the power-law nature of p/cr^, one could ask 
whether or not the (3 vs. 7 behavior is independent of scale-free pjo'^ . Since there is no known 
theoretical motivation for either of these relationships, we attack this question with different tactics. 

The extended secondary infall model (ESIM) (Williams, Babul, &: Dalcanton 2004; Austin 
et al. 2005) is a semi-analytical halo formation scenario based on the work of Ryden &; Gunn 
(1987). In ESIM, spherical halos form by accreting shells of material that have decoupled from 
Hubble expansion. These halos exist at overdense locations and their formation includes effects 
from secondary perturbations that can be completely specified and controlled. As shells collapse, 
the shells' energies change in response to the continually changing potential. In fact, the ESIM 
halo collapse process is best described as violent relaxation because direct two-body effects do not 
take place (there are no shell-shell interactions), while the shells are allowed to exchange energy by 
interaction with the global potential. In general, ESIM halos have power-law density distributions 
over the radial ranges in which N-body simulations show decidedly non-power law behavior, e.g., 
where NEW profiles change from p oc to p oc r~^. 

Standard ESIM halos have scale-free p/o"^ (Austin et al. 2005; Barnes et al. 2005), but what 
do their P vs. 7 curves look like? Panel a of Eigure 11 represents the P and 7 values taken 
from a standard ESIM halo, and the locus bears no resemblence to the linear relation found in 
N-body simulations. This implies that the P-j relationship is not just a manifestation of a scale- 
free p/<J^. It is certainly possible that the process responsible for creating scale- free p/cr'^ is not 
fully expressed in standard ESIM simulations, leading to a different /3~7 relationship. However, 
we suggest that the near-linear relation between /3 and 7 supports an earlier argument that one 
generic outcome of mostly-radial collapses is an instability that produces velocity distributions 
with nearly isotropic cores and more radially anisotropic envelopes (Barnes et al. 2005). Due to 
the semi- analytical nature of the ESIM formalism, this instability cannot develop in ESIM halos. 
However, by manipulating the secondary perturbations present in ESIM halo formation, one can 
induce this type of anisotropy profile in ESIM halos. Doing so radically changes the associated 
density profile from the single power-law of the standard halo to an NFW-like form (Barnes et al. 
2005). These changes strongly affect the overall shape of the P vs. 7 plot. Panel b of Eigure 11 
illustrates the P vs. 7 relationship for an NFW-like ESIM halo. While the points certainly do not 
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fall on a line, they more closely approximate the N-body simulation results. The implication is 
that the physics driving the density profile shape is also important to establishing the linear (3-^ 
relationship. 

6. Summary & Conclusions 

It is well-established empirically that p/a'^ is a power-law in radius for virialized, collisionless 
halos (Taylor & Navarro 2001; Austin et al. 2005; Dehnen & McLaughlin 2005). We exploit this 
relationship and transform the usual Jeans equation into the constrained Jeans equation, solely in 
terms of the logarithmic density slope 7 and the anisotropy (3. Doing this allows us to solve directly 
for 7 once an anisotropy distribution has been chosen. 

In an earlier paper, Barnes et al. (2006), we discussed the solutions of this constrained Jeans 
equation when the anisotropy is assumed to be zero everywhere. Here, we have extended our earlier 
work to include constant, but nonzero, anisotropies as well as anisotropy distributions. In particular, 
we have looked at two flexible, but distinct, types of anisotropy distributions; a power-law and a 
hyperbolic tangent. We find that, in general, the solutions of the constrained Jeans equation that 
most resemble the empirical NFW and N04 density profiles have [3\o ~ and /3hi ^ 0.5. Anisotropy 
must be an increasing function of radius in halos like the ones produced in N-body simulations. 

A central conclusion of our present work is that there is a distinct similarity between the 
shapes of radial density and anisotropy profiles. Density profiles that best match the NFW form 
are found when the anisotropy profile is of the tanh form (right panels of Figure 7) , which contains 
the Osipkov-Merritt and Mamon-Lokas forms as special cases. N04-like solutions are found with 
either profile (right panels of Figures 4 and 8), but with the tanh [3 profile, the anisotropy must be 
roughly constant over the range 10~^rvir < r < r^\^ (Figure 8). 

Another aspect of this profile similarity is given by the (3 vs. 7 plots (Figures 5, 6, 9, and 10). 
The combination of (3 and 7 profiles that produce the most linear curves (as suggested by the results 
of N-body simulations; Hansen &: Moore 2004, 2006), are, of course, those profiles that resemble 
each other most closely. The inflection point in the tanh anisotropy distribution also appears in the 
NFW density slope distribution (right panels of Figure 7 and Figure 9). In other words, when this 
second-derivative symmetry is lost, as with NFW- like density and power-law anisotropy profiles 
(right panels of Figure 3), so is the near-linear correlation between (3 and 7. However, we note that 
the most linear (3-^ relations found here are also those that are the most quantitatively different 
from the results of simulations, in that the solid red lines in Figures 6a and 9a extend outside the 
simulation uncertainty range. This suggests that the relationship between (3 and 7 for any single 
halo is not strictly linear. 

Utilizing the results of semi-analytic models of halo formation, we argue that the (3-^ relation- 
ship is not just a manifestation of a scale-free p/cj'^ ■ Rather, a near-linear relation between [3 and 
7 supports an earlier argument regarding the nature of halo density profiles (Barnes et al. 2005). 
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Velocity distributions with nearly isotropic cores and more radially anisotropic envelopes result 
from radial collapses. This form of anistropy distribution necessarily creates a density profile that 
is less steeply rising in the central regions, since there are fewer radial orbits that plunge close to the 
center and can support a strong cusp. Likewise, more tangential orbits in the outer regions would 
increase the density since the average time an orbit spends far from the center would increase. At 
least in a qualitative sense, this picture supports the observed similarities between profile shapes; 
more isotropy (smaller f3) implies less steeply varying density (smaller 7) and vice versa. 

The outcome of this work is that there appear to be at least two physical processes at work 
in the formation of collisionless dark matter halos in N-body simulations. The evidence presented 
suggests that one process creates the scale- free /o/cr^ distribution and must be a fairly generic 
collapse process as the p/cr^ behavior is seen in the more restrictive ESIM halos. As a consequence, 
another process links the anisotropy and density profiles in a unique, nearly linear way. 
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Fig. 1. — An overview of the models and results presented in this paper. The columns refer to the 
two types of density profiles discussed here; Navarro-Frenk- White (Navarro, Frenk, & White 1996, 
1997, NFW) and Navarro et al. (2004, N04). These profiles are described in detail in §2. The rows 
delineate the three types of anisotropy profiles; constant (§3), power-law (§4.1), and hyperbolic 
tangent (§4.2). This table is intended to guide the reader through and summarize the various 
figures and key results. 
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Fig. 2. — The 7 = —d\n{p/pQ)/d\n{r/rQ) distributions (solid red lines) that result from solving 
the constrained Jeans equation (Equation 7) with constant anisotropy (3 values compared with 
N04 (Equation 6, dash-dotted green lines) and NEW (Equation 5, dashed blue lines) 7 profiles. 
Panels a-c show the solutions that best-fit the N04 profile; no combinations of (3 and a produced 
solutions that approximate NEW profiles. Panel a shows the solution with a = 1.875, (3 = 0.1. The 
solution in panel b has a = 35/18,/? = —0.1, and panel c has a = 1.975, /? = —0.2. The solution 
in panel d has a = 1.875, /3 = 0.5. It is included to illustrate the behavior of radially anisotropic 
systems. All solutions, regardless of a values, with (3 > 0.3 resemble this solution and provide poor 
approximations to any empirical 7 distribution. See §3 for details. 
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Fig. 3. — With power-law I3{x), panels a, c, and e show the calculated 7 profiles (solid red lines) that 
best match the NFW 7 profile (dashed blue lines) for a = 1.875, 35/18, and 1.975, respectively. 
To magnify the variations, each 7 profile has been normalized by the solution. The dot-dashed 
green lines show the N04 7 profile. Panels b, d, and f illustrate the anisotropy distributions used 
to calculate the solutions (solid red lines). The orange triple-dot-dashed lines are scaled versions 
of the solution 7 profile and are included to compare the shapes of the P and 7 profiles. The thick 
dotted lines represent an anisotropy distribution commonly used to describe the results of N-body 
simulations (Mamon & Lokas 2005, Equation 60). The vertical dotted lines mark the boundaries 
of the fitting region, x = 0.01 and x = 10. The solutions never well-approximate the NFW profile. 
In fact, for a = 1.875, the solution is much more like an N04 profile. Note that an NFW 7 profile 
has an inflection point at Inx = 0, whereas this /3 proflle does not. See §4.1 for details. 
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Fig. 4. — With power-law (3{x), panels a, c, and e show the calculated 7 profiles that best match 
the N04 7 profile for a = 1.875, 35/18, and 1.975, respectively. Panels b, d, and f illustrate the 
anisotropy distributions used to calculate the solutions. The line styles are the same as in Figure 3. 
The solutions match the designated N04 profile reasonably well, in contrast with the results shown 
in Figure 3. 
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Fig. 5. — Correlations between power-law velocity anisotropy /? and logarithmic density slope 
—7 profiles of constrained Jeans equation solutions for 10~^2;vir < x < Xvh- Panels a, b, and c 
represent the best-fit NFW solutions with a = 1.875, 35/18, and 1.975, respectively. The dashed 
blue line represents the mean trend found in simulation results (Hansen &: Moore 2006) and the 
dash-dotted green lines mark the extent of variations from those simulations. While the solutions 
are quantitatively consistent with the simulations, their shapes (particularly the concavity) are 
unlike the results of simulations over similar ranges of 7 values. 
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Fig. 6. — Correlations between power-law (3 and —7 solutions for 10~^Xvir ^ x < Xvir- Panels 
a, b, and c represent the best-fit N04 solutions with a = 1.875, 35/18, and 1.975, respectively. 
The dashed blue and dash-dotted green lines are the same as in Figure 5. The most linear (and 
hence most qualitatively consistent with simulations) relationship in panel a is inconsistent with 
the quantitative results of simulations. As in Figure 5, the red solution curves that quantitatively 
agree with simulations have more concavity than the simulations themselves. 
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Fig. 7. — With tanh panels a, c, and e show the constrained Jeans equation solution 7 profiles 
that best match the NFW 7 profile for a = 1.875, 35/18, and 1.975, respectively. Panels b, d, and 
f illustrate the anisotropy distributions used to calculate the solutions. The line styles are the same 
as in Figure 3. Unlike the results shown in Figure 3, NFW profiles are reproduced very well by these 
solutions. Note that now the 7 and /? profiles share the inflection point behavior at Inx = 0. This 
similarity of 7 and /? profile shapes implies a near linear correlation in the /3-7 plane (Figure 9). 
See §4.2 for details. 
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Fig. 8. — With tanh /3(x), panels a, c, and e show the calculated 7 profiles that best match the N04 
7 profile for a = 1.875, 35/18, and 1.975, respectively. Panels b, d, and f illustrate the anisotropy 
distributions used to calculate the solutions. The line styles are the same as in Figure 3. These 
solutions are fair approximations of the designated N04 profile. However, note that the /3 profiles 
are nearly horizontal. Higher a values lead to moderate variations in anisotropy, but not as large 
as those in Figure 7. 
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Fig. 9. — Correlations between tanh /3 and —7 for 10~^Xvir ^ x < Xvh- Panels a, b, and c represent 
the best-fit NFW solutions with a = 1.875, 35/18, and 1.975, respectively. The dashed blue and 
dash-dotted green lines are the same as in Figure 5. As in Figure 6, the most linear relation is 
quantitatively inconsistent with the results of simulations. Higher a values produce red curves that 
have larger regions of linearity and are consistent with simulations. Recall that in this case the 7 
and (3 profiles have similar shapes (particularly the inflection point at Inx = 0). 
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Fig. 10. — Correlations between tanh /3 and —7 for 10~^Xvir < 2; < Xvir- Panels a, b, and c represent 
the best-fit N04 solutions with a = 1.875, 35/18, and 1.975, respectively. The dashed blue and 
dash-dotted green lines are the same as in Figure 5. Overall, the red curves are more horizontal 
than the simulation results. However, larger a values allow for quantitative agreement between the 
solution curves and simulations. 
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Fig. 11. — The (3 vs. —7 curve for the standard ESIM halo (a). The corresponding curve for the 
ESIM halo that has been modified so that the density profile resembles that of an NFW model. 
The points in each panel marked by red solid circles represent parts of the halo that lie in the 
range 10~^ < r/r2oo < 0, where r2oo is the virial radius of those simulations. The open circles lie 
beyond the virial radius. The dashed blue and dash-dotted green are the same N-body simulation 
results shown in the previous plots of (3 vs. 7. Despite the bump at 7 = 2, panel b resembles the 
results from N-body simulations much more closely than does panel a. Each of these ESIM halos 
has power-law p/cj'^profiles, but they have radically different 7-/? relationships. This suggests that 
the mechanisms that give rise to each of these relations are independent. The standard halo (panel 
a) is nearly isotropic out to the virial radius while the NFW-like halo (panel b) has a tanh-like 
P profile. This prompts us to postulate that the P-j connection is rooted in the instability that 
creates isotropic cores and radially anisotropic outer regions, and that the p/a^ relationship is due 
to a more general physical process. 



